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Abstract; In the present work, a dynamic model is 
developed for the Multi-effect evaporator (MEE) to 
study the transient behavior of the system. Each 
effect in the process is represented by some variables 
which are related to the energy and material balance 
equations for the feed, product and liquor flow. 
Backward feed is used for the development of the 
model for six effect evaporator system. For the steady 
state and dynamic simulation, the ‘fsolye’ and 
‘ODE4S’ solvers in MATLAB source code is used 
respectively. 


Keywords: Multi-Effect Evaporator, backward feed, 
fsolve, ode45. 


L. INTRODUCTION 


Some researchers have solved the mathematical models 
of MEE system of various process industries by using 
different solution techniques to obtain steady state 
simulation [3-20], Some of them also have performed 
steady state simulation for MEE system of the paper 
industry [1], [2], [3]. [9-11], [16-17], [20] & [29-34] by 
solving the models using numerical techniques, The 
dynamic behavior of two or three effect MEE systems in 
process industries like sugar, food, desalination and 
paper etc. is studied by few researchers [13] and [21-26]. 
Stefanov and Hoo [27] developed a distributed 
parameter model of black liquor falling film evaporator, 
based on first principles knowledge of fluid dynamics 
and heat transfer processes for evaporation side of the 
single plate (lamella) of a falling film evaporator. 
Stefanov and Hoo [28] further expanded a single plate 
model to develop a fundamental model of a falling film 
evaporator that accounted for the condensation side of 
the plate, the heating/flashing at the evaporator entrance, 
the evaporator inventory, mixing and recirculation flows 
but neglecting the effect of BPR. 


In the present investigation, an attempt has been made for 
the study of dynamic responses for tubular type falling 
film MEE system. The lumped parameter model (System 
of ordinary differential equations) of sextuple effect 
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falling film evaporator system of a paper industry is 
developed by using unsteady state energy and material 
balance equations and Physico-thermal properties of 
black liquor including the effect of boiling point rise for 
backward feed sequence, The transient behavior of the 
system is studied by disturbing the input parameters like 
feed flow rate, feed concentration, feed temperature and 
steam temperature. For the steady state and dynamic 
simulation ‘fSolve’ and ‘ode45’ solver in MATLAB 
source code is used respectively, 


Il, MATHEMATICAL MODELLING 


The Mathematical modelling is carried out for sextuple 
effect falling film evaporator system with backward feed 
In backward feed evaporator system, the steam input is 
given in the first effect and the feed liquor input is in the 
last effect as shown in the Fig. 1. 

Model equations are developed for ith effect using 
material and energy balance equations ¢ 
The equations stating the physical properties of the black 
liquor are taken from [32]. It is assumed that the vapor 
generated by the process of concentration of black liquor 
is saturated, It is also assumed that the energy and mass 


shown in [32]. 


accumulation in the vapor lump is very small as compared 
to the enthalpy of the steam and can be neglected. 


Material balance for liquor in the i" effect: 


uu MI, (t) = WI,,, — WI, — Wy, 
dt 


Energy balance for liquor in the i" effect: (2) 


(x4 (hl )=WI, hl, -WI,hl, —Wy,hy, + Wy, hv, -Wy. he 
It 


where 


where Qis the rate of heat transfer and also 
Q=U,A,(TV,.,-Tv, -BPR,) 


shell area of the evaporator and Li is the liquor level fori” 


Material balance for solids in ith effect: 
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a imi 
The vapor and liquor in ith effect are in equilibrium and 
the relation for the liquor and vapor temperature is 
defined in terms of boiling point rise (BPR) is as follows: 
T1,=Tv, + BPR, 4) 


where boiling point rise (BPR) is a function of 
temperature and solid concentration 

Differentiating equation (4) with respect to time we get 
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MI, can be written as: 
MI, = AL,p, (6) 


where pi is the density of the liquor for ith effect, which is 
the function of temperature and solid concentration, A is 
the effect. Differentiating equation (6) with respect to 
time 
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Fig. 1: Flow diagram of Sextuple Backward feed evaporator 


Since pi is a function of temperature TI and concentration 

X, this equation reduces to equation (8) by substituting of 
4.5 4) from equation (5) and rearranging the terms. 
dt 
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Comparing the resultant differential equation (8) with 
equation (1) we get 
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c4 = aL 4{ 2 \ Zone )+(% 
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Enthalpy is a function of temperature Differentiating 
d 
=~TI(t) 
MA()*hL(1) with respect to time and substituting dt 
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from equation (5) and dt from the equation (8) 
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Comparing the resultant differential equation (10) with 
equation (2) we get 
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Finally differentiating M1(t)*X(t) with respect to time 
and substituting the value of a M_, (t) from equation (8) 
and rearranging the equation and representing it in the 
form of coefficients 


where 
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Thus the equations (9), (11) and (12) form a set of 
differential equations representing the liquor flow inside 
the evaporator and the solid content of the liquor for the 
ith effect. These final equations, represehts the dynamics 
of multi-effect evaporators, derived from basic energy 
and mass balance equations and using various terms 
related to physical properties of black liquor, water and 
steam like density and specific heat of the liquor, boiling 
point rise (BPR), enthalpy of water and saturated steam. 
Further equations are applied to all the six effects for i= 
1,2,3,4,5and6, 


II. STEADY STATE SIMULATION 


For obtaining the steady state solution of the sextuple 
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effect evaporator system there is a need to solve the 
system of nonlinear simultaneous equations given by (9), 
(11) and (12) under steady state conditions, The detailed 
study of steady-state simulation is given in our earlier 
work Kumar etal. [31], 
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IV. DYNAMIC SIMULATION 


For the dynamic simulation first order nonlinear 
differential equations (8), (9) and (10) are solved 
simultaneously in the same order for all the six effects of 
the backward feed sextuple effect evaporator system. 
Steady-state solution of the model provides the initial 
values of the system variables at time, t= 0. The Solution 
of such types of simultaneous nonlinear ordinary 
differential equations is extremely intricate in nature even 
by using sophisticated numerical techniques. In the 
present investigation, an attempt has been made for 
dynamic simulation for tubular type falling film MEE 
system with split feed sequence by using “ode45’ solver in 
MATLAB source code respectively. 


V. RESULTS AND DISCUSSION 
A, Effect of varying feed flow rate 


Since MEE system is in backward feed, hence the input 
liquor enters the evaporator system from the last effect. 
Thus for the variation in the feed flow rate the disturbance 
is applied in the last effect. The effect of +10% step input 
in feed flow rate on the temperature and concentration of 
last (6") and first (1") effects are shown in the Fig. 3 to Fig. 
4 respectively. The temperature and concentration of both 
the effects show an increase or decrease with decrease or 
increase in the feed flow rate. It is obvious as fresh steam 
supply rate is constant and water to evaporate and 
decreases per unit time. The Time constant (TC) for each 
plot is obtained and shown in the plot. TC indicates that 
response time for temperature change is much less than 
that of concentration change, Also, the steady state in 6th 
effect is achieved at a faster rate than the Ist effect. This is 
noticeable as in backward feed sequence; feed enters first 
in the 6th effect, so the disturbance dampened out in the 
6th effect at the faster rate, 


Variation in the Temperature of 6th effect with 10% decrease in feed flow rate 


60 7 T 
55 
= 50 
i 
45 
40 = 
a 500 1000 15002000 
Time (Sec) 


\Variation in the Temperature of 6th effect with 10% increase in feed fow rate 
60 7 T , ; = 


Variation in the Product Concentration of {st eect with 10% decrease in feed flow rate 


iy 
© 50h j 
F a 
5} 4 
40. 
0 500 Yoo = 1500» 200025003000 
Time (Sec) 
Time (Sec) 


\Vatiation in the Product Concentration of 6th effect with 10% decrease In fead flow rate 
0.16 T 


04 i: f nl atl 
05 1 15 2 25 3 35 4 

Time (Seo) eat 
\Vatiation in the Product Concentration of 6th effect with 10% inerease in feed flow rate 
0.18 7 t T T ———T 


2 
= 


A warKa) 


Tim Bex) nd 


Figure 3 Response of 6th effect by disturbing + 10% in the 
feed flow rate 
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Figure 4. Response of Ist effect by disturbing + 10% 
in the feed flow rate 


B. Effect of varying feed concentration 


For the variation in the feed concentration, the disturbance 
is applied in the last effect. The effect of +10% step input 
in feed concentration on the temperature and 
concentration of last and the first effect are shown in the 
Fig. 5 to Fig. 6 respectively. The dynamic behavior of 
effect’s temperature on disturbances in feed concentration 
shows a slight but insignificant change in temperature. 
However the change. it observed is unidirectional i.e. the 
temperature increases irrespective of increase or decrease 
in feed concentration. The changes in product 
concentration of both the effect show increase or decrease 
according to as the feed concentration is increase or 
decrease. This is because the fact that AT across the 
evaporator system remains constant and vapor-liquid 
equilibrium of each effect remains almost unchanged for 
the optimum performance. Time constant shown in the 
figures indicates that response time for temperature 
change is much less than that of concentration change. 
Since the feed is backward thus, the steady state in 6th 
effect is achieved ata faster rate than the Ist effect. 
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Figure 5. Response of 6th effect by disturbing + 10% 
in the feed concentration 
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Figure 6. Response of Ist effect by disturbing + 10% in the 
feed concentration 


C. Effect of varying steam temperature 


In backward feed MEE system, live steam enters in the 
first effect. Thus for the variation in the steam 
temperature. the disturbance is applied in the first effect. 
The effect of £10% step input in steam temperature on the 
temperature and concentration of last and first effects are 
shown in the Fig. 7 to Fig. 8 respectively. The 10% change 
in steam temperature results in increase or decrease in the 
temperature of both the effects before obtaining the steady 
state for 10% increase or decrease respectively. 10% 
disturbance in steam temperature does not result in any 
noticeable change in the product concentration of both the 
effects. However after scaling down Y-axis value, it was 
observed that product concentration increase and then 
decrease or decrease and then increase for 10% increase or 
decrease in the steam temperature respectively, 
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Figure7. Response of 6th effect by disturbing + 10% 
in the steam temperature 


‘Variation in the Temperature of 1st etfect with 10% decrease in steam temperature: 
190 


ox Te = 306 sec 


D. Effect of varying feed temperature 


Since feed liquor enters in the last effect, thus for the 
dynamic response of the feed temperature the disturbance 
is applied in the last effect. The effect of +10% step change 
in feed temperature on the temperature and concentration 
of last and the first effect are shown in the Fig. 9 to Fig. 10 
respectively. It is evident from the figures that 10% 
disturbance in feed temperature does not bring noticeable 
change in the temperature and product concentration each 
effect. However after scale down Y-axis, it is observed 
that temperature of the both the effect increases and 
decreases to obtain the steady-state with an increase and 
decrease in feed temperature and the product 
concentration of each effect first decreases and then 
increases to obtain the steady state with a very small 
fluctuations about the steady-state up to four to five 
decimal places in the value of concentrations of both the 
effect and conversely for 10% increase in feed 
temperature. 


\Verietion in the Temperature of 6th ettect with 10% decreane in feed temperature 


° 500 1000 4500 


Time (Seo) 


72000 2500 


Variation in the Temperature of 6th effect with 10% Increase in feed temperature: 
a “r 


ey g 
8 110 g 
= 
400 
“ 0 an) soo —-1000~=~=«NSOO.~=~COS=S~*«RSCSC OOO 
° a a a a 5 te 
Time (Seo) Variation in the Produot Concentration of eth etfest with 10%» decrease in feed temperaturt 
Variation |n) the Temperature of 1st elfect with 10% increase in steam temperature 016, v “ : 
130, 
Bor | 
120. 3 
g = 12 
& 110f rs 
f on aa 
oo) TC =390 see 0 500 ~~1000~=~*~«SSD~S~COSCOSCSC OO 
sol ‘ | Time (See) 
° ooo ——«1500«=« 000500 K 
Variation in the Conaentration of 6th affect with 10% increase in feed temperature 
ae o.16¢-— dibs ai : 


Variation in the Concentration of 1st effect with 10% decrease in steam temperature 


X (kg’kg) 


Te = 21107 800 I] 
oo 
6 7 8 o 
x10" 
Variation in the Product Concentration of 1st effec! with 10% increase in steam temperature 


45 
Time (Sec) 


~ — 
[T= em a | 
oOo. 2 s «¢ 6 6 7 8 
‘Time (Seo) 
Figure 8. Response of Ist effect by 


disturbing + 10% in the steam temperature 


=a 500 000 1500 


Time (Sec) 
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Figure 10. Response of Ist effect by disturbing 
+ 10% in the feed temperature 


VI. CONCLUSION 

Dynamic modeling is an useful tool in the determination 
of process variables in transient conditions. For 
designing an effective and efficient control system, it is 
always desirable to have a thorough knowledge of the 
system behavior in different conditions as well as 
detailed knowledge of variation in the process variables 
under transient condition. There is always a de 
process industry for easier to use and more affordable 
advanced control technologies and products. Dynamic 
model helps us in better understanding of the process and 
behavior of its variables, thereby helps in determining a 
better control system. Thus in the present investigation, 
an unsteady-state lumped parameter models were 
developed for sextuple effect falling film evaporator for 
backward feed sequence for concentrating the black 
liquor by using material, energy balance equations and 
parametric correlations. For the steady state and dynamic 
simulation, the *fsolve’ and *ode45’ solvers.in MATLAB 
lidated for the steady 
state solution using literature data. The dynamic 
behavior of each effect’s temperature and product 
concentration was studied by disturbing the liquor flow 
rate, feed concentration, steam and feed temperatures by 
=10%, The transient study shows that the steady state is 
reached more quickly for temperature in comparison of 
the solid concentration and all of the responses converge 
inasmooth fashion. 
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